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Abstract 

It has been shown by Strahov and Fyodorov that averages of products and ratios of cha- 
racteristic polynomials corresponding to Hermitian matrices of a unitary ensemble, involve 
kernels related to orthogonal polynomials and their Cauchy transforms. We will show that, 
for the unitary ensemble 4- | det M \ 2a e - nV( - M ^dM of nxri Hermitian matrices, these kernels 
have universal behavior at the origin of the spectrum, as n — > oo, in terms of Bessel functions. 
Our approach is based on the characterization of orthogonal polynomials together with their 
Cauchy transforms via a matrix Riemann-Hilbert problem, due to Fokas, Its and Kitaev, and 
on an application of the Deift/Zhou steepest descent method for matrix Riemann-Hilbert 
problems to obtain the asymptotic behavior of the Riemann-Hilbert problem. 

1 Introduction 

Characteristic polynomials of random matrices are useful to make predictions about moments of 
the Riemann-Zeta function, see [5l ll8llT^ll21j . Another domain where they are of great value is 
quantum chromo dynamics, see for example [21 El E3 E3j • In this paper we consider characteristic 
polynomials det (a; — M) of random matrices taken from the following unitary ensemble of n x n 
Hermitian matrices M, cf. 15} I25j 

J-|detM| 2a e- ntrV/(M) dM, a>-l/2. (1.1) 

Here dM is the associated flat Lebesgue measure on the space of n x n Hermitian matrices, 
and Z n is a normalization constant. The confining potential V in Ijl.lj) is a real valued function 
with enough increase at infinity, for example a polynomial of even degree with positive leading 
coefficient. This unitary ensemble induces a probability density function on the n eigenvalues 
x\ , . . . , x n of M, see j^H] 



1 n 

P (n) (xi, ...,x n ) = — Y[w n (xj)A 2 (xi, . . . ,x n ), 



where A(xi, . . . ,x n ) = Yli<j( x j ~ x i) stands for the Vandermonde determinant, where Z n is a 
normalization constant (the partition function), and where w n is the following varying weight 
on the real line 

Wn {x) = \x\ 2a e- nV ( x \ (1.2) 
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The unitary ensemble is relevant in three-dimensional quantum chromo dynamics 33 , 
and has been investigated before in [H 151 12U1 1251 128| . where universal behavior for local eigenvalue 
correlations is established in various regimes of the spectrum, as n — > oo. 

It is known that averages of products and ratios of characteristic polynomials are intimately 
related to orthogonal polynomials and their Cauchy transforms, see |H1 El Ell 130] , Let 
7Tj tn (x) = x 3 + • • • be the j-th. degree monic orthogonal polynomial with respect to w n . There is 
an integral representation for the monic orthogonal polynomials, which appears already in the 
work of Heine in 1878, see for example 31 , 



7Tn,n( X ) = / • ■ ■ / Y\( X ~ X (xx, . . .,X n )dxi . . . dx r . 
J J ,=1 



So, the monic orthogonal polynomial 7r„ in can be understood as the average of the characteristic 
polynomial det(x — M) over the unitary ensemble (jl.ip 

(det (a; - M)) M = K n>n {x). 

Here, the brackets are used to denote the average over the ensemble (|1.1|) of random matrices M. 
A first generalization of this formula was obtained by Brezin and Hikami jS] , and also by Mehta 
and Normand |27| . They have derived a determinantal formula for the average of products of 
characteristic polynomials in terms of orthogonal polynomials. A further generalization was 
obtained by Fyodorov and Strahov |17j , who derived a determinantal formula for the average of 
both products and ratios of characteristic polynomials in terms of both orthogonal polynomials 
and their Cauchy transforms. Here, the ratios gave rise to the Cauchy transforms. For explicit 
formulas and streamlined proofs of these results we refer to (7j. 

Recently, Strahov and Fyodorov (HOI showed, see also for an alternative proof, that the 
averages of characteristic polynomials of n x n Hermitian matrices, are governed by kernels 
related to orthogonal polynomials and their Cauchy transforms 

hj n (z) = — !— f 7Fj ' n X w n (x)dx, for Imz / 0. (1.3) 
2tti J x — z 

Namely, kernels Wi in+m made of orthogonal polynomials, kernels Wn tn+m made of both orthogo- 
nal polynomials and their Cauchy transforms, and kernels Win n+m made of Cauchy transforms 
of orthogonal polynomials. See Tabled for the explicit expressions of these kernels. This con- 
nection between the averages of characteristic polynomials and the three kernels is given by, see 

[a ism 

(ndet(^-M)det( K -M)\ = A * ^ det (Tfr,M*fo,yj)) 1<iJ<t , 

\i=l / M ilj=n C i,n \ > 



n 



det(y; — M) \ 



\ det(xi — M) j 



h0t=±l f , k A(x,y) 



(-1) 2 ( Cn _l jn ) 



M 



A 2 (x 



y^27^y det ( w ii,*( x i>yj))i<ij<h 



and 



n det(Xj _ M )/ M 



(-iy 



k \pn—k—X,n) 



\k n-1 



n 



det (W/jj, n _fc(a; <r ( i ),x .( fc+J -))) 



l<i,j<k 



( 2k )- ^ n _ k a f^ &{Xo(l),----,Xo{k))^{x (j{k+1) ,...,X (j{ 2k)y 
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Finite kernels 


W I)n+m (C,r]) 


1"n+m.n(C)l"n + m-l,n(»)) — 71"n + m-l,n(C) 7r n+m.n('?) 


C-v 


W n ,n+rn((,V) 


h„ + m,n(0 7T n + m-l,n(v) — hn + m-l,n (0 7r n + m,n (v) 


C-v 


WlII,n+m(G,lj) 


^■ra + m.n (O^ri + m.— l,n ^n + m — 1 ,n (O^n + Tt.n ( if) 


C-ri 



Table 1: Expressions for the finite kernels Wj jn+m ,Wji jn+m and Win,n+m, cf. |HU] . 



where x = (xi, . . . , x^), y = (y\, . . . , y^), where Cj )U = —2mjj n with jj >n the leading coefficient 
of the j-th degree orthonormal polynomial with respect to w n , and where S2k is the permutation 
group of the index set {1, . . . , 2k}. There are also explicit formulas for averages containing non- 
equal number of characteristic polynomials in the numerator and the denominator, in terms 
of these kernels, see |3f)j for details. Strahov and Fyodorov pHl used this connection, together 
with the Riemann-Hilbert (RH) approach, to establish universal behavior, as n — > oo, for the 
averages of characteristic polynomials of random matrices taken from the unitary ensemble 

J_ e -ntrV(M) dMj (L4) 

in the bulk of the spectrum. 

It is the goal of this paper to establish universal behavior as n — > oo, for the kernels 
Wj jn+m ,Wji jn+m and Wju jn+m (and thus also for the averages of characteristic polynomials) 
associated to the unitary ensemble 1)1. ljl . appropriate scaled at the origin such that the asymp- 
totic eigenvalue density at the origin is 1. This scaling limit is called the origin of the spectrum 
by various authors, see for example )HIHl l2L)l I25j . It will turn out that this universal behavior is 
described in terms of the Bessel kernels given in Table 01 For the case a = 0, our results agree 
with those of Strahov and Fyodorov [SO] . 

The issue of universality at the origin of the spectrum for the averages of characteristic poly- 
nomials, corresponding to Hermitian matrices of the unitary ensemble 1)1.1)1 . was also considered 
by Akemann and Fyodorov They showed, on a physical level of rigor using Shohat's method, 
that the asymptotic behavior near the origin, as n — > oo, of the orthogonal polynomials and 
their Cauchy transforms are expressed in terms of Bessel and Hankel functions, see |Bj for details. 
However, explicit expressions for the universal behavior of the three kernels Wj jn+m ,Wji jn+m 
and Wni,n+m at the origin of the spectrum have not been given yet, which we will determine 
on a mathematical level of rigor using the RH approach, as in [M)\ . 

In jS] was assumed that the potential V is an even polynomial with positive leading coeffi- 
cient, and that the spectrum support is only one interval. In this paper, we can allow V to be 
quite arbitrary, and assume the following 

V : R -> R is real analytic, (1.5) 
V(x) 

lim - — z ^ • — = +oo, (1.6) 

|a;|-»oc log(x 2 + 1) 

^(0) > 0, (1.7) 
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Limiting Bessel kernels 



a+ i a+ i J a+ iW a -i(^)-J a -iK)./ a+ i(^) 
^ 71 W^) 



1 iH-J? (1) iW)ViN 



2 W=v) 



(2) 



n 11 s(C-v) 



-7rC 2 r? 



8(C-»?) 



^ 77 8[C^ 



Case a = 



sin 7r(f — 77) 
i"(C-»?) 



j e i *"(C-'?) 

2T(C-r;) 



j e -i7r(C-i?) 

" 27r(C-r?) 







Table 2: Expressions for the limiting Bessel kernels. Here, J v is the usual J-Bessel function of 
order u, and and are the Hankel functions of order v of the first and second kind, 
respectively. The right column denotes the expressions in case a = 0. 



where ip is the density of the equilibrium measure fj,y in the presence of the external field V, 
see ^JEHUHl- The equilibrium measure /ly has compact support, it is supported on a finite 
union of intervals (since V is real analytic), and it is absolutely continuous with respect to the 
Lebesgue measure, i.e. dfiy(x) = ip(x)dx. The importance of fj,y lies in the fact that its density 
tp is the limiting mean eigenvalue density of the unitary ensemble . Assumption (|1.7|) then 
states that the mean eigenvalue density does not vanish at the origin. 

Our results are given by the following three theorems. We use C+ and C_ to denote the 
upper and lower half-plane, respectively. 

Theorem 1.1 Fix m G Z, let Wi^ n + m be the kernel given in TableUi and let 7,- n > be the 
leading coefficient of the j-th degree orthonormal polynomial with respect to w n . For £,77 G C 

7 ™ +m - 1 '"^(o) Wl ' n+m G^o) ' ^0)) 

= (m/>(0)) 2a e ny ^ (e%M {C+v) S a ,l{(,ri) +0(l/n)J , (1.8) 

as n — > oo, where the Bessel kernel S a j(C, rj) is given in Tabled The error term holds uniformly 
for C and rj in compact subsets of C. 

Theorem 1.2 Fix m G Z, Zet W/j. n + m 6e £/ie kernel given in Table an<i 7j jft > be 
the leading coefficient of the j-th degree orthonormal polynomial with respect to w n . Then the 
following holds. 
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(a) For C £ C+ and ijfC 



2 ^ _r7 w f C 

7n+m-l,n~jJ7^\ W II,n+m 



= (C-^)e"^ (C "" ) JI+ // (C^)+0(l/n), (1.9) 

as n — ► oo ; where the Bessel kernel //(Ci 7 ?) 75 given in Tabled The error terra holds 
uniformly for £ and 7/ in compact subsets o/C+ and C, respectively. 

(b) For C € C_ andrj£C 

2 (-V ixr ( C g 

T«+«-i.» n ^(0) ^^+™ ^( ) ' n ^(0) 

= (C - ^e-^T^-^J- J7 (C, r?) + 0(1/7-0, (1-10) 

as n — > oo ; where the Bessel kernel 3Ztt(C>v) 7S given in Tabled The error term holds 
uniformly for £ and n in compact subsets of C_ and C, respectively. 

Theorem 1.3 Fix m G Z, /ei Win. n+m be the kernel given in Table QJ and let 7j,n 

> 6e 

t/ie leading coefficient of the j-th degree orthonormal polynomial with respect to w n . Then the 
following holds. 

(a) For (,r/G C+ 

7 ' +m - 1 ' n ^(0) WlII > n+m (^0) ' ^0) ) 

(1 \ ^ a / v'(o) \ 

^yj e- nV{0) [e-^ {C+V) ^,ni(C,v)+0(l/n)j , (1.11) 

as n — ► oo, where the Bessel kernel -^///(Ci 7 /) 7S given in Tabled The error term holds 
uniformly for £ and r/ in compact subsets o/C+. 

f&j For C G C + and n G C_ 

2 1 rxr Z' C " 

Vm-l,n~T7^ W /n,n+m 



' nV'(O) JJJ >" +m ^"^(0) ' n^(0) 



i \ 2a 



nv(o) f e -S§}«+")j± m (C,7/) +O(1/t0) , 



;i.i2) 



as n — > oo, where the Bessel kernel ^ miCiV) 75 given in Tabled The error term holds 
uniformly for £ and 77 in compact subsets of C + and C_ , respectively. 

(c) For C,77 e C_ 

7 " +m - 1 '™^(0) ^ m ' n+m (^0) ' ^0) ) 

(1 \ ^ a f v'(o) \ 

^yj ^^i (C+7?) J- m (C^) + 0(l/n)J , (1.13) 

as n — ► 00, where the Bessel kernel ///(Ci 7 /) 7S given in Tabled The error term holds 
uniformly for £ and 7/ in compact subsets of C_ . 
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Remark 1.4 In case a = we can simplify the expressions for the limiting Bessel kernels, using 
the facts that, see ^ 



Ji(z) = \ — sinz, J_i(z) = \ — cosz, H^\z) = —i\ — e iz , 
2 V vrz 2 V vrz 2 \ nz 

H {1 1 (z) = i/^e* Hf\z) = iJ^e^, (z) = \f^e~^. 

~2 V TTZ 2 V TTZ ~2 V VTZ 

We then obtain the kernels given in the right column of Table [21 This is in agreement with the 
results of Strahov and Fyodorov [HO]- Note however that in [HO] the second and the third kernel 
are multiplied with an extra factor —27ri. 

Remark 1.5 As noted before, it has been shown by Strahov and Fyodorov [30], see also [7], 
that 



2m ^hm WlI ' n G 



\ de t(-C_ J -M)/ M '""^^(O) JJ ' n V^(0)'n^(0); • 

Then it follows from ()1.9j) . Tableland ^ formula 9.1.3], that for ( G C+, 



de , _ / 2 ^ + |K)^-|K)-^-iK)^ + |K))+od/-)> 



as n — > oo, where Y"^ is the Bessel function of the second kind of order v. By formula 9.1.16], 
the right hand side of this equation is 1 + 0(l/n), as it should be. Similarly we find the same 
result for £ G C_. 

The proofs of these theorems are based on the characterization of orthogonal polynomials 
with respect to the weight Q1.2|) . together with their Cauchy transforms via a 2 x 2 matrix 
RH problem for Y, due to Fokas, Its and Kitaev ^H], and on an application of the Deift/Zhou 
steepest descent method ^3] for matrix RH problems. See |1UI I22j for an excellent exposi- 
tion. This technique was used before by Deift et al. ^Qj to establish universality for the local 
eigenvalue correlations in unitary random matrix ensembles (jl.4j) in the bulk of the spectrum. 
Strahov and Fyodorov j20j used this method also to establish universality for the three kernels 
W^n+m, Wji jn+m and Wjn jn+m in the bulk of the spectrum. 

In a previous paper [^3J together with A.B.J. Kuijlaars, the asymptotic analysis of the RH 
problem for Y, corresponding to the weight (|1.2|) . has already been done. An essential step in 
the analysis is the construction of the parametrix near the origin, which gives us the behavior of 
Y near the origin. Here, the Bessel functions come in. In [23], the behavior of the first column 
of Y (with the orthogonal polynomials as entries) was determined near the origin for positive 
(real) values, and used to establish universality for the local eigenvalue correlations at the origin 
of the spectrum, in terms of a Bessel kernel. Here, we determine the behavior of the first column 
of Y, as well as the second column of Y (with the Cauchy transforms of orthogonal polynomials 
as entries) in a full neighborhood of the origin, and use this in a similar fashion to prove our 
results. 

The rest of the paper is organized as follows. In Section [2] we give a short overview of the 
asymptotic analysis of the corresponding RH problem for Y . In Section [H] we determine the 
behavior of Y near the origin, in terms of Bessel functions. This will be used in the last section 
to prove our results. 
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2 The corresponding RH problem 



In this section we recall the matrix RH problem for Y, due to Fokas, Its and Kitaev |16j . which 
characterizes the orthogonal polynomials with respect to the weight (|1.2j) . together with their 
Cauchy transforms. We also give a short overview of the Deift/Zhou steepest descent method 
to obtain the asymptotic behavior of Y. For details we refer to ^3 EH], see also [TfHHI]. 

Our point of interest lies in the asymptotic behavior, as n — > oo, of the orthogonal polyno- 
mials itn+m,n of degree n + m with respect to the weight w n , for any fixed m £ Z. So, in contrast 
to the RH problem in |13l I25| , we have to modify the asymptotic condition at infinity of the RH 
problem, and leave the jump condition unchanged. However, this will not create any problems. 
We seek a 2 x 2 matrix valued function Y = y( n + m > n ) that satisfies the following RH problem, 

cf. [ig E3 EH OSES]. 

RH problem for Y: 

(a) Y : C \ R -» C 2x2 is analytic. 

(b) Y possesses continuous boundary values for x £ R \ {0} denoted by Y + (x) and Y-(x), 
where Y + {x) and Y— (x) denote the limiting values of Y(z') as z' approaches x from above 
and below, respectively, and 

Y+(x) = Y„(x)r m \ , forxGl\{0}. (2.1) 

(c) Y has the following asymptotic behavior at infinity: 

Y(z) = (I + 0(l/z))\^ Q z _ (n+m) J, asz^oo. (2.2) 

(d) Y has the following behavior near the origin: 

'1 \z\ 2a " 



Y(z) 



O ( 1 j z j 2a ) , if a < 0, 
0{] J') . iCn >U. 



(2.3) 



as z -» 0, z £ C \ R. 

Remark 2.1 The O-terms in condition (d) of the RH problem are to be taken entrywise. So for 

example Y(z) = O I , , 2q , ) means that Fu(js) = O(l), Fi 2 (z) = 0(|z| 2q ), etc. This condition 

is used to control the behavior of Y near the origin. In the following we will not go into detail 
about this condition, and refer to |23l I32j for details. 

The unique solution of the RH problem for Y, see |16| (for condition (d) we refer to |23j). is 
then given by 

/ ^"n+m. n {%") ^n+m,ra(^) \ 

Y(z) = Y { - n+m ' n \z) = , (2.4) 

\~^™ln+m-l,n ,T ri+m-l,n\Z) —2' n "i'ln+rn-l,nh'n+m-l,n{z) J 

for z £ C\K, where 7Tj in is the j-th degree monic orthogonal polynomial with respect to w n , 
where jj n is the leading coefficient of the j-th degree orthonormal polynomial with respect to 
w n , and where h^ n is the Cauchy transform of 
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Remark 2.2 The superscript n + m in the notation y( n + m -™) corresponds to the asymptotic 
condition (c) at infinity of the RH problem, which yields that the orthogonal polynomials in 
the solution (|2.4jl of the RH problem have degree n + m and n + m — 1. The superscript n 
corresponds to the jump condition (b), which yields that the orthogonality is with respect to 
w n . 

Remark 2.3 We note that the first column of Y contains the orthogonal polynomials, and the 
second column their Cauchy transforms. So, from Tableland (|2.4j) . the kernel Wi tn+m depends 
only on the first column of Y, the kernel W//,n+m on both the first and the second column, and 
the kernel Wni,n+m only on the second column, as follows: 



7n+m-l,n -27Ti(C 



and 



W '"^ K - "» = 7 ; + l,„-^K-» det Gwo ' (2 - 7) 

The asymptotic analysis of the RH problem for Y includes a series of transformations Y i— > 
T i — > S i — > to obtain a RH problem for i? normalized at infinity (i.e. R(z) — > / as n — > oo), 
and with jumps uniformly close to the identity matrix, as n — > oo. Then |1U1 1131 HI] , i? is also 
uniformly close to the identity matrix, and by unfolding the series of transformations we obtain 
the asymptotic behavior of Y. 

Before we can give an overview of the series of transformations, we need some properties of 
the equilibrium measure [iy for V. Here, we closely follow [23], see also j!21 113j . The support 
of fiy consists of a finite union of intervals, say \ < Jj=x~[bj-i,aj], and we define its interior as 
J = \Jj^i~(bj—i, dj). The N + 1 intervals of J are referred to as the bands. The density ip of the 
equilibrium measure is given by 

^( x ) = J-R^ 2 ( x )h( x ), for x G J, (2.8) 

with h real analytic on R, and where R is the 2(./V + l)-th degree monic polynomial with the 
endpoints a^, bj of J as zeros, 

JV+l 

R(z)= Hiz-b^z-aj). (2.9) 

i=i 

We use i? 1 / 2 to denote the branch of V^R which behaves like z N+1 as z — ► oo, and which is 
defined and analytic on C \ J. In l|2.8|) . ii^ 2 is used to denote the boundary value of R 1 / 2 on J 
from above. The equilibrium measure satisfies the Euler-Lagrange variational conditions, which 
state that there exists a constant I G K such that 

\og\x - s\ip(s)ds-V(x) =£, for x & J, (2.10) 

2 /" log |z - s\i(j(s)ds - V(x) < £, forxGM\J. (2.11) 

If the inequality in (|2.11|) is strict for every i£l\J, and if h(x) ^ for every x G J, then V is 
called regular. Otherwise, there are a finite number of points, called singular points of V, such 
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that h vanishes there, i.e. a singular point in J, or such that we obtain equality in (|2.11j) . i.e. a 
singular point in R \ J. 



Let (T3 = ( o -i ) be the Pauli matrix. Following JH]) see also [25] . we define the 2x2 matrix 
valued function 



T(z) = e -(™+ m )| CT 3y( z ) e ("+ m )| CT 3 e -(™+m)3Wa 3 1 for z £ C\ 



(2.12) 



where I is the constant that appears in the Euler-Lagrange variational conditions (12.10(1 and 
([2.11]) . and where the scalar function g is defined by 



g(z) = / log(z — s)ip(s)ds, for z£C \ (— oo, ajv+i]- 



(2.13) 



Note the small difference in the definition ([2.12]) of T with its definition in [1.3] 125] , which comes 
from the modified asymptotic condition (c) of the RH problem for Y. For the case m = 0, both 
definitions agree. It is known [13] 125] that T is normalized at infinity and satisfies the jump 
relation T + (x) = T^{x)v^ l \x) for x G R \ {0}, where 



v {1) (x) 



' e -(n+m)(g + (x)-g-(x)) \ x \2a e mV(x) 

e (n+m)(g+(x)-g-(x)) 



x E J\{0} 



-2Tri(n+m)Qj I \2a ^mV(x) g(n+m)(g+ (x)+g_ (a;)— V(a;)— N 







, x G (aj,bj), 



'\ \x\ 2a e mV ^ e^ n+rn ^ g+ ^ +9 -^~ v 
.0 1 



, x < bo or x > &jv+i • 



(2.14) 



The constant J7j is the total /iy-mass of the iV + 1 — j largest bands. 

The second transformation is referred to as the opening of the lens. Define [25] for every 
z G C \ M lying in the region of analyticity of h the scalar function 



(*) 



^ r a N+i 



R 1/2 (s)h(s)ds, 



(2.15) 



where the path of integration does not cross the real axis. Then |25j . on the bands, <j) is purely 
imaginary and satisfies 



20+ (x) = -2<p_(x) = g + (x) - g-(x), 



for x £ J, 



(2.16) 



so that 2(f> and — 2<p provide analytic extensions of g+ — g- into the upper half-plane and lower 
half-plane, respectively. The opening of the lens is based on the factorization of the jump matrix 
i/ 1 ) on the bands, see (|2.14j) . into the following product of three matrices, cf. |25j 



-(n+m)(g + (x)-g_(x)) 



\ x \2a e mV(x) N 
,(n+m)(g+ (x)-g-{x)) 



1 N 

\ x \-2a e -mV(x) e -2{n+m)4>- (x) -y 



\x\ 2a e mV ^\ ( 1 s 

\ x \-2a e -mV(x) q J I I ^,1 ~2a e ~mV (x) e -2{n+m)4> + (x) ^ 
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Figure 1: Part of the contour S. 



We take an analytic continuation of the factor |x| 2Q e ml/ ( :r ' 1 by defining for z in the region of 
analyticity of V, 



(-z) 2a e mV( - z \ 



if Re z < 0, 
if Re z > 0, 



(2.17) 



with principal branches of powers. We now open the lens. Let S be the lens shaped contour, 
as shown in Figure E going through the endpoints aj, bj of J, going trough the origin, and also 
going through the singular points of V in J. Of course we take the lens shaped regions to lie 
within the region of analyticity of 4> and V. 



Define, cf. 



S(z) 



T(z) 



T(z) 



T(z) 



- W ( z )-l e -2(n+m)<^) 1 



oj(z)~ 1 e~ 2( - n+m ^ ( - z ^ 1 



for z outside the lens, 

for z in the upper parts of the lens, ^ ig\ 
for z in the lower parts of the lens. 



As for the first transformation F h T, there is small difference in the definition IJ2.18J) for S 
with its definition in |25| . which comes from the modified asymptotic condition (c) of the RH 
problem for Y. For the case m = 0, again both definitions agree. Then [2Sj, the matrix valued 
function S is normalized at infinity and satisfies the jump relation S+(z) = S-(z)v^ (z) for 
z£S, where 



v^(z) 



1 



t o \ z \2a> e mV{z) 



, ze£nc±, 



-2a -mV{z) 







, zeJ\{0}, 



' e -2ni(n+m)£lj i^i 2a e mV(z) e (n+m)(g+(z)+g- (z)-V(z)- 
n p 2ni(n+m)£lj 



(2.19) 



, z€(a,j,bj) 



\ \ z \2a e mV(z) e (n+ m )(g + (z)+g-(z)-V(z)-^ 



,0 



1 



, z < bo or z > ajv+i- 



For z in a neighborhood of a regular point x G J we have, cf. |25| . 
Re 0(z) > 0, if Im z / 0, 
and for every regular point in R \ J we have from the Euler-Lagrange variational condition 

CUD, cf. [HI 



5+(x) + g-{x) - V(x) - I < 0, for x G R \ J. 
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So, we expect that the leading order asymptotics are determined by a RH problem for p(°°) , nor- 
malized at infinity, that satisfies the jump relation P^°\x) = P^°\x)v^°°\x) for x 6 (6o,ajv+l)j 
where 



v {oo) (x) 



q | 2 ,|2o e m,y(x) 
^ x ^-2a e —mV(x) q 



, forxGJ\{0}, 

(2.20) 

for x € (dj, bj), j = 1 . . . N. 



The solution of this RH problem is referred to as the parametrix for the outside region, and is 
constructed using a Szego function on multiple intervals associated to \x\ 2a e mV ( x \ cf. [2*5] . and 
using Riemann-Theta functions, cf. JB], see also For our purpose here, we do not need the 
explicit formulas for p(°°), and refer to |13l I25j for details. 

Before we can do the third transformation, we have to be careful since the jump matrices 
for S and p(°°) are not uniformly close to each other near 0, near the endpoints a,i,bj of J, 
and near the singular points of V. To solve this problem, we surround these points by small 
non-overlapping disks, say of radius 5 > 0, and within each disk we construct a parametrix P 
satisfying the following local RH problem. 

RH Problem for P near xq where xq is 0, an endpoint of J, or a singular point of V: 

(a) P(z) is defined and analytic for z £ {\z — xq\ < 5q} \ X for some 5q > 5. 

(b) P satisfies the same jump relations as S does on S n {\z — xq\ < 5}. 

(c) There is k > such that, as n — > oo, 

P(z) (p(°°A 1 (z)=I + 0(l/n K ), uniformly for \z - x \ = 6. (2.21) 



(d) SP 1 has a removable singularity at xo- 

For regular endpoints and the origin we can take k = 1 in (|2.21|) . It is known that this local 
RH problem is solvable for every xq. For the endpoints of J and the singular points of V we 
refer to for the origin we refer to |25| . For our purpose here, it suffices to know the explicit 
formula for the parametrix near the origin. 

We will now give the explicit formula for the parametrix P near the origin, see |251 Section 
5] for details, see also |321 Section 4]. This is an essential step in the asymptotic analysis of the 
RH problem since it allows us to determine the behavior of Y near the origin, which will be the 
main tool to prove our results. Introduce the scalar function 

f(j\ _ / ^(s) -*^+(0)» iflmz>0, . . 

J{Z) \ -uf>{z)-i<f> + (Q), iflmz<0, 1 > 

which is defined and analytic in a neighborhood of the origin. The behavior of / near the origin 
25 Section 5] is given by 

f(z) = irip(0)z + O(z 2 ), asz^O. (2.23) 

Let Us be the disk with radius 5 around the origin, with 5 > sufficiently small such that Us lies 
in the region of analyticity of 4> and V. Since /'(0) = irip(0) > we can choose 5 also sufficiently 
small such that / is a conformal mapping on Us onto a convex neighborhood of 0. We have that 
fix) is real and positive (negative) for x £ Us positive (negative). 
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Figure 2: Decomposition of Us and f(Us) into eight regions. 



Decompose f(U$) into eight regions I- VIII, as shown in the right of Figure El divided by 
eight straight rays 



7T . 



Tj = {C E C | argC = (j - j=l v . . ,8. 

This in turn divides the disk Us into eight regions I'-VIIF as the pre-images under / of I— VIII, 
as shown in the left of Figure El Sector I' and IV correspond to the right and left upper part of 
the lens inside Us, respectively, sector V and VIII' to the left and right lower part of the lens 
inside Us, respectively. 

Let ^ a be the piecewise analytic matrix valued function (221 Section 4], see also [23 Section 
5], that satisfies the jump relation ^ a>+ (() = \P a) _(£)u a (£) for ( £ U^?, where 



Va(C) 



1 
-1 



1 

— 2ma i 

1 

-,2-nia j 



for CeTiUTg, 

for c g r 2 u r 6 , 
for c g r 3 u r 7 , 
for c g r 4 u r 8 , 



and that has the following behavior near the origin, 

icr icr 



if a < 0, and 

*«(C) = 



icr icr 

icr ici- 
icr icr 



o 



as C —> 0, 



asC^Ofor | < |argC| < ^, 



ICI" 

icr 



icr 



, as C — > for < | arg (\ < f and ^ < | arg C| < vr 



if a > 0. The behavior of *$> a near the origin will ensure that part (d) of the RH problem for P 
is satisfied, see (23 for details. The matrix valued function ty a is constructed out of Bessel 
functions of order a ± | , and its explicit formula for < arg £ < f is given by 



*«(0 = ^C 1/2 



(2.24) 
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b H a } bj a J+l b J+x 

Figure 3: Part of the contour The points z\ and zi are singular points of V . 
For j < arg (, < ^ it is given by 

*" (C) = rrVi T tr ~k ,1,2 J I -* e ^™' (2 ' 25) 

where 1^ and ifj, are the modified Bessel functions of order v. See j521 Section 4] for the explicit 
expressions of f Q in the other sectors of the complex plane. Also define the piecewise analytic 
function W by 

f z a e mY r 1 , if z g iir,iv',v',vr, 

W(z) = (2.26) 
[ (-z) Q e m — , if G I' ,11', VII', VIIF. 

And finally, define the following matrix valued function, analytic in a neighborhood of the disk 

U s , 



E nW (z) = E (z)e^ n+m ^^e-^^= ft ! ) , (2.27) 



where the matrix valued function E is given by [251 (5.27)-(5.30)]. 

Then, cf. |25| Section 5], the parametrix near the origin is defined by 

P{z) = E n+m Az)^ a ((n + m)f(z))W(z)-^e-^ +m ^ z ^. (2.28) 

Remark 2.4 In contrast to |25| Section 5], we evaluate the matrix valued function *$> a in 
(n + m)f(z) instead of in nf(z). This comes from the fact that, in order that the matching 
condition (c) of the RH problem for P is satisfied, we need to cancel out the factor e -( n + m )H z ) (7 3 
instead of e~ n ^ z \ This follows in essence from the modified asymptotic condition (c) of the RH 
problem for Y. For the case m = 0, the definition (|2,28|) of the parametrix P near the origin 
agrees with its definition in |25| Section 5]. 

Now, we have all the ingredients to give the third transformation. Define |13l I25j the 2x2 
matrix valued function R as 

f S(z)(P^°°yy 1 (z), for z outside the disks, 

R(z) = < ' (2.29) 

[ S(z)P~ 1 (z), for z inside the disks. 

Then |13l I25j . R is normalized at infinity, and analytic on the entire plane except for jumps 
on the reduced system of contours £r, as shown in Figure 01 and except for possible isolated 
singularities at the endpoints aj, bj of J, at the singularities of V and at 0. However, from 
condition (d) of the RH problem for P, these singularities are removable, so that R is analytic 
on C \ Sj?. It is known ^3 [25] that the jumps of R on are uniformly close to the identity 
matrix as n — > 00. This implies ^Hl> see also |10| 114) 

R(z) = I + 0(l/n K ), asn->oo, (2.30) 

uniformly for z G C \ where K is the constant that appears in the matching condition (c) of 
the RH problem for P. By tracing back the steps Y \— > T S ^ R we obtain the asymptotic 
behavior of Y in all regions of the complex plane, as n —* 00. 
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3 Behavior of Y near the origin 

In this section we unravel, as in |25l Lemma 7.1], the series of transformations Y i— > T S i— > R, 
see Section to determine the behavior of the first and the second column of Y inside the disk 
Us- This behavior will be the main tool to prove our results. Note that the second column of Y 
has jumps on the real axis, see (|2.1jl . So, for the behavior of the second column of Y inside the 
disk Us we have to distinguish between the upper and lower parts of Ug. 

For notational convenience we introduce the 2x2 matrix valued function, cf. |25l Lemma 

7.1] 

M(z) = M n+mjn (z) = R(z)E n+m , n (z), for z G Us, (3.1) 

where E n+rrijn is given by 1)2. 27j) . For the case m = 0, the M-matrix defined by 1)3.1)1 corresponds 
to the M-matrix in |25l Lemma 7.1]. It is known that M is analytic on Us, that each entry of 
M is uniformly bounded in Us as n — > oo, and that detM = 1, cf. |251 Lemma 7.1]. 
We also need the following lemma. 

Lemma 3.1 For z G Us, 

2g(z) -2<P(z) -£= V(z). (3.2) 

Proof. Let H(z) = 2g(z) — 2<p(z) — £ — V(z), which is defined and analytic for z G Us\M.. For 
x £ (-5, 5) C J we have by lf2~ToT> 

H + (x) = H_(x)=g+(x)+g„(x)-£-V(x), (3.3) 

so that H is analytic in the entire disk Us- For x G (—8,5) we have by (|2,13|) 

g+(x) + g~(x) = 2 / log |s — s| , 0(s)cis. 



Inserting this into (|3.3|) and using the Euler-Lagrange variational condition ()2.1fl|) , we have that 
H(x) = for x G (—5,(5). This implies from the uniqueness principle that H = on L^, which 
proves the lemma. □ 

First, the behavior of the first column of Y inside the disk Us is given by the following 
theorem. 

Theorem 3.2 Fix m G Z. For z G Us and n sufficiently large, the first column ofY = y( n + m > n ) 
is given by 

(^J) = z^e^V^e-^e^^Miz) 

'((n + m)f(z)) 1/2 J a+h ((n + m)f(z)) 
^n + m)f(z)) 1/2 J a _ 1 _((n + m)f(z)) 

Here, J u is the J -Bessel function of order v, f is given by (|2.22|) . and M is given by ()3.1j) . 



(3.4) 



Proof. Let z be in sector I' of the disk Us, see Figure|2j Unfolding the series of transformations 
Y^T^S^Rwe obtain by (l2~T2]) . (|^TH|) . (|2~27?|) and (l2~2lH) 



Y(z) = e( n+m ^^R(z)E n+m:n (z)y a ((n + m)f(z))W(z) 



-<T3 



x ^(n+m)^ f ( ^ rle i (n+m) ^ (2) Jje-(^)^ e (^)^. (3.5) 
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Note that u(z) = z 2a e mV{z \ see (t2~TTl) . and that W(z) = {-z) a e mY ^ = ^ e -™ e mZ ^ ) see 
(|2.2fij) . Inserting this into (|3.5|) and using (|3.1|) and (|3.2|) , the first column of Y is then given by 

= z ~ aen ^ e{n+m) ^ 3M ( z )^<*(( n + m )f( z )) e7TiaC73 (J) • ( 3 - 6 ) 

Since f(z) is in sector I of f(U$), see Figured we have for n sufficiently large (namely n + m > 0) 
that < arg(ra + m)f(z) < 7r/4. So, we have to use (|2.24|) to evaluate W a ((n + m)f(z)). From 
(|3.6|) and 1, formulas 9.1.3 and 9.1.4], which connect the Hankel functions of the first and second 
kind with the ordinary J-Bessel functions, we then establish 1)3 .4(1 in sector I' of Us- 

Now, let z be in sector IF of Us- Similarly as in sector F, we obtain by (|2.12|) . (|2.18|) . (|2.28j) 



and (I2T29|) 

Y{z) = e^^R(z)E n+m , n (z)^ a ((n + m)f(z)) 

X W ~(z)~ as e -( n + m )<t > ( z ) a 'A e ~(n+m)^a 3 e (n+m)g(z)ai _ 

Since W(z) = z ^ e -^ e m^ ^ gee and using (j^jj) and the first column of Y is then 

given by 

= z- a e nZ ¥e {n+m ^ a3 M(z)* a ((n + m)f(z))e mar73 Q . (3.7) 

Since ir/4 < arg(n + m)f(z) < ir/2 for n sufficiently large, we have to use 1)2.25(1 to evaluate 
^a{{n + m )f{z))- This implies, using pfl formula 9.6.3], which connects the modified Bessel 
function I a ±i with the Bessel function J a ±i, that 



V a ((n + m)f(z)y ia ° 3 



_ _ ((n + m)f(z)) 1/2 I a+ i((n + m)f{z)e ") 

7re 2 | i/2 „ | (3.8) 

i((n + m)f(z)) I a _i({n + m)f(z)e~-) 



7TJ. 

\/ire * 



'((n + m)/(z)) 1/2 J a+l _ ((n + m)/(z)) X 
((n + m)/(z)) 1/2 J Q _i((n + m)/(z)) 



2 

Inserting this into (|3.7|) . we establish ()3.4I) also in sector IF of Us- 

In the other sectors of the disk Us the calculations are similar, and are left as an easy exercise 
for the careful reader. □ 

Remark 3.3 For the case m = 0, this theorem agrees with |25l Femma 7.1]. 

Remark 3.4 It is not quite clear from (|3,4j) that the first column of Y is analytic in the entire 
disk Us, which must be the case since it has polynomials as entries, see 1)2.4)1 . Obviously, it is 
analytic on Us \ (—6, 0]. From 1, formula 9.1.35] we have for x G (—5, 0) 



1 /2 "* 

-(-(n + m)f(x)) J a+ i(-(n + m)f(x)) 
(_(„ + m)/(x)) 1/2 J Q _ i (-(n + m)/(x)) 



(3.9) 



So, the first column of Y is analytic in the entire disk Us except for a possible isolated singularity 
at the origin. Since J a +±(z) = 0(z a=t 2) as z — > 0, see |TJ formula 9.1.10] this singularity is 
removable, which implies that the first column of Y is analytic in the entire disk. 
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Next, the behavior of the second column of Y in the upper part of the disk Us is given by 
the following theorem. 

Theorem 3.5 Fix m G Z. For z £ (/jfl C+ and n sufficiently large, the second column of 
Y = y( ra + m > ra ) i s given by 

= ,2 e n 2 -0re 4 e l + ; 2 3 M(2r) 
\*22(*)/ 2 



({{n + m)f(z)) 1/2 H®i ((n + m)/(z)) N 
((n + m)/(z)) 1/2 i/« , ((n + m)/(z)) y 



(3.10) 



i?ere ; is the Hankel function of the first kind of order u, f is given by (|2.22j) . and M is 

given by (|3.1|) . 

Proof. Let z be in sector I' of Us- Instead of (|3.6j) we get for the second column of Y 

= zQe " n ^ e(n+m) ^ 3M ( z )^(( n + m )/( z )) e7riQ<J3 (J) • ( 3 - n ) 



Since < arg(n + m)f{z) < 7r/4 for n sufficiently large, we have to use (|2.24j) to evaluate 
^a({ n + m )f( z ))- Inserting this into (|3.11j) . we obtain (|3.1U|) for this choice of z. 

Now, let z be in sector II' of Us- Instead of (|3.7|) the second column of Y is given by 

^ 12 jjj) = z a e- n ^e ( - n+m ^ a3 M{z)y a ((n + m)f(z))e 7Tiaa:i (f \ . (3.12) 

Since ir/4 < arg(n + m)f(z) < ir/2 for n sufficiently large, we have to use (|2.25p to evaluate 
^a\{ n + m )f( z ))- From ^ formula 9.6.4], which connects the modified Bessel function K a± i 

with the Hankel function H, , of the first kind, we then have 

«.4 



1 Trig 

e 2 

7T 



1 



1 ((n + m)f(z)) 1/2 K a+h ((n + m)f(z)e-f) 
:i((n + m)f(z)) 1/2 K a _i ((n + m)/(*)e~^r) 

'((n + m)/(^)) 1/2 lf« i ((n + ro)/(*))' 
v ((n + m)/(z)) 1/2 i^ ((n + m)/(^ 



Inserting this into ()3.12j) . equation (|3.1U|) is proven in this sector as well. 

Similarly, we can prove ([3.10)1 in the other sectors of the upper part of Us- □ 

And finally, the behavior of the second column of Y in the lower part of the disk Us is given 
by the following theorem. 
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Theorem 3.6 Fix m G Z. For z G Us D C_ and n sufficiently large, the second column of 
Y = y( n + m > n ) j s given by 

= - z "e- n ^^e-^e( n+m ^M(z) 




({{n + m)f{z)f l2 H^ +1 _ ((n + m)/(z)) X 
((n + m)f{z)f l2 H^\ ((n + m)/(z)) 



(3.13) 



Here, Hi, is the Hankel function of the second kind of order v, f is given by (|2.22[) . and M is 
given 6y (|3.1|) . 

Proof. The proof is similar to the proofs of Theorem 13.21 and Theorem 13.51 □ 
Remark 3.7 By Q2.1|) . the jump relation for the second column of Y is 

" few) ^ (Sw) »^ n ' K fo " eR ^ (3 ' 14) 

For x G (0,5) one can check easily, using ()3.2j) . (|3.1Ujl . ()3.13|) and [IJ formulas 9.1.3 and 9.1.4] 
that (ETT^I is satisfied. For x G (-5,0) it follows from (pOSjl . (I3~TT)|) . (I3T31 and 1, formulas 9.1.3, 
9.1.4 and 9.1.39] that (l3~Tl|) is satisfied. 

Remark 3.8 Theorems 13.21 15*31 and l3~^l give because of ()2.4|) . after straightforward calculations, 
the behavior near the origin of the orthogonal polynomials and their Cauchy transforms. It has 
been shown before by Akemann and Fyodorov 6 that the behavior near the origin of the 
orthogonal polynomials is given in terms of the J-Bessel functions and that the behavior 

near the origin of their Cauchy transforms is given in terms of the Hankel functions H f_ 1 of the 

(2) 

first kind in the upper half-plane, and in terms of the Hankel functions H ! 1 of the second kind 

in the lower half-plane. However, in [B] this was done on a physical level of rigor, and under the 
assumption that the eigenvalue density was supported on only one interval. 



4 Proof of Theorem fiTTI LTTTTI 

In this section we prove the universal behavior at the origin of the spectrum for the three 
kernels Wj jn+m , Wn,n+m and Wni,n+m, i n terms of the Bessel kernels given by Table[2J Similar 
as in |24| I25j . where we have investigated local eigenvalue correlations, we do this by using the 
connection of these kernels with the solution of the RH problem for Y, see (|2.5|l — (|2.7|) . and by 
using the behavior of Y near the origin, derived in the previous section. 

We first need the following lemma's. 

Lemma 4.1 Let M be the matrix valued function given by (|3.1D . and let (,?/£C. Then, each 
entry Mij of M satisfies 

uniformly for ( and rj in compact subsets ofC. 
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Proof. Let (,f/ £ C, denote Cn = ^jm and = w ^ - ) , and let 7 be a positively oriented 
simple closed contour in Us going around the origin. Then, since M is analytic on Us, an 
application of Cauchy's formula shows that 

Mij(z) 



MijiCn) ~ Mijirjn) = (Cn - Vn)^: j 



(z - Cn)0 - Vn) 



dz, 



for £ and n in compact subsets of C and n sufficiently large. Since M{j is uniformly bounded in 
Us as n — > 00, see Section |3 the integral is uniformly bounded for £ and rj in compact subsets 
of C as n — > 00. This proves the lemma. □ 

Lemma 4.2 Fix m G Z. Le£ £,n G C ; and denote 

Cn = (n + m)f ( -— ) , and fj n =(n + m)f 



T/ien, 

r°Cy 2 «/ Q± x(C„) = r Q (vrC) 1/2 J a±| (7rC) + 0(l/n), (4.2) 

as n — > 00, uniformly for £ in compact subsets of C. TTie /e/£ /land side 0/ (|4.2|) is uniformly 
bounded for C in compact subsets of C as n — > 00. ^4/so 



= r a K) 1/2 «/ a± iK) - "-"(vrn) 1 / 2 ^!^) + ( 
as n — » 00, uniformly for £ and n in compact subsets of C 



(-rj 



n 



(4.3) 



Proof. By (f2~2^ it follows that C„ = vrC(l + 0(l/n)) as n -> 00, uniformly for £ in compact 
subsets of C. Inserting this into the left hand side of ()4.2|) we easily obtain estimate (|4.2(l . 
cf. [23 Lemma 7.2]. Since d^(C) = ("Hv(0 with ^ entire, see p] formula 9.1.10], we have 
that C _a+5 d n ,,i (0 is entire. This implies by (14.21) that the left hand side of (14.21) is uniformly 
bounded for £ in compact subsets of C as n — > 00. 

Let K\,K2 be compact subsets of C, and let 7 be a positively oriented simple closed contour 
with K\ and K2 in its interior. Define 

q n (z) = z~ a ~z l J 2 J a±] _{z n ) - z- a {Ttz) l / 2 J a±} _(Kz), (4.4) 

with z n = (n + m)/( - ). Note that g n is analytic in an open neighborhood of the interior of 
7 for n sufficiently large. An application of Cauchy's theorem then shows that 

UC) - qn (v) = (C - 17)^ (z-CK^-.) ^ 

for £ G -fTi and n G and n sufficiently large. Since q n {z) = 0(l/n) as n — > 00 uniformly for 
z in compact subsets of C, see (|4.2j) and (|4.4I) . and since £ and n are not close to the contour 7, 
the lemma is then proven. □ 

Now, we are ready to prove the universal behavior at the origin of the spectrum for the 
kernel Wj n+m . 
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Proof of Theorem II. 1L Let (,r] G C, denote 



Cn 



mp(0) ' 



Vn 



, Cn = {n + m)f(Cn), and fj n = (n + 



and let Y = y( n + m > n ). Similar considerations as in |24[ I25j . using (|2.5j) and the behavior ()3.4|) 
of the first column of Y inside the disk Us, show that, 

W hn+m (C,T]) = 7n+m-l,n L W/, ra+m (Cn, %) 



nV(0) 

All (Cn) ^11 fa*) 
27rz(C-r?) V*2l(Cn) *2lfan) 



det 



(nV(0)) 2a e?( y ^) +v '^)) 



1 



x det 



M(Cn 



c- a Cn /2 j„ +i (c 



C a Cn /2 J n _l(( ri 



+ M{r} r , 



The matrix in the determinant can be written as, cf. |241 I25j . 



M(Cn) 



r a Cn /2 J Q _l(Cn) 




7r a rjn /2 J a _i(rj, 



(4.5) 



+ M(C„)- 1 (M( ??n )-M(C„)) 



r/~ 
r/- 



^1/2 T 
'/n ^o-I 




Since det M = 1 and each entry of M is uniformly bounded in Us as n — » oo, see Sectional each 
entry of M(£ n ) _1 is uniformly bounded for £ in compact subsets of C as n — > oo. Using Lemma 

1 /2 

14.11 and the fact that i]~ a fj n J a ±i(fj n ) is uniformly bounded for 77 in compact subsets of C as 



n 



00, see Lemma 14.21 we then obtain 



M(C„)- 1 (M(r/ n )-M(C„)) 



7] a fjl /2 J a+ i(Vn) 
rf -^ 2 J a _i{fj n ) 



~l/2 

Using the facts that det M = 1 and that C~ a Cn J a ±i((n) is uniformly bounded for £ in compact 
subsets of C as n — > 00, see Lemma 14.21 we then find 

W Iin+m (C,v) = (n^(0)) 2a e^ v ^+ v ^ 




1 



det 



C a Cn /2 J a+ . 



(Cr, 

(C, 



.1/2 



V Vn « Q+ . 



2 



+ 0(l/n) 



(4.6) 



We can now replace z a Zn 2 J + \(z n ) by z a {iTz) l l 2 J , i(7rz) for 2 = C,r/, and obtain the 

uzn 2 win 2 

limiting Bessel kernel $ a ,i(CiV) given in Table El However, then we make an error which does 
not hold uniformly for £ and r/ close to each other. To solve this problem we will work as in 
|24| I25j . We subtract the second column in the determinant from the first one. Prom (|4.H|) and 
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1/2 

the fact that i]~ a ri n J a ±i(fj n ) is uniformly bounded for 77 in compact subsets of C as n — ► oo, 
the term inside the brackets in Q4.6J1 is then given by 

i ( r a «) 1/2 ^ + |K) - rw% } w ry-^y 2 j a+l (^ n )\ 

Using the fact that 

^n /2 ^ ± l(^)=^ Q (^) 1/2 ^±l(^)+0(Vn), 

and the fact that 

C - 77 

remains bounded for £ and 77 in compact subsets of C, which follows since z~ a (irz) 1 ^ 2 J , 1 (z) 
is entire, we then easily obtain that the term inside the brackets in ()4.6|) is given by 



det :_ Q : u/2 r* j ' « ' 1/2 r 5 ' +o(v»). (4.7) 



2(C-r?) lr Q K) V2j a-iK) ^"W^-lM 



2 2 



The first term in (|4.7|) is exactly the limiting Bessel kernel Ja,/(C)77), see Table |2 From (j4.6j) 
and (j4.7|) we then obtain 

W Ijn+m (C,v) = (nmf a e% {V{(: " )+VM) (Sa,i((,v) + 0(l/n)y (4.8) 
as n —* oo, uniformly for £ and 77 in bounded subsets of C. Note that 

e ?(V(f„)+V(^)) = e ^(0)+^(C+r?)+O(l/n) 

= e ny We^T (C+??) (l +0(l/n)), asn^co, 

uniformly for £ and 77 in compact subsets of C. Inserting this into 1)4. 8j) . and using the fact that 
^ajiCiV) is bounded for £ and 77 in compact subsets of C the theorem is then proven. □ 

In order to prove Theorem 11.21 we also need the following lemma, which is analogous to 
Lemma 14.21 

Lemma 4.3 Fix m £ Z. Let £ € C+, and denote Q n = (n + m )f( ni p(o) )- Then, 

e^H^n) = C(0 1/2 H^(7rC)+0(l/n), (4.9) 

as n ^ 00, uniformly for (" in compact subsets o/C+. TTie Ze/t /land side 0/ ()4.9|) is uniformly 
bounded for £ in compact subsets of C+ as n ^ 00. 

Proof. Recall, cf. the proof of Lemma 14.21 that Q n = (ttQ(1 + 0(l/n)). Inserting this into 
the left hand side of 1)4.9)1 and using the fact that is analytic on C \ (— oo,0] we easily 

obtain estimate (|4.9jl . Since is analytic on C \ (-00, 0], we have that CiO 1 ^^^ «) 

is bounded for £ in compact subsets of C\ (— 00, 0], and thus in particular in compact subsets of 
C+. Together with estimate (|4.9[) this implies that the left hand side of 1)4.9)1 remains uniformly 
bounded for £ in compact subsets of C+ as n — > 00. □ 



20 



Proof of Theorem OJ Let ( G C+, r? G C, denote 



Cn 



n^(0) ' 



tin 



, Cn = (n + m)/(Cn), and r/ n = (n + m)f(r] n ), 



and let Y = y( n + m > n ). Instead of equation Q4.5JI we obtain from (|2,6|) . from the behavior (|3.4|) 
of the first column of Y inside U$, and from the behavior ()3.1(Jj) of the second column of Y in 
the upper part of Us, 

Wn,n+m(C,V) = ln+m-X,n ~ , W H ,n+m(Cn, Vn) 



nip(0) 

(Y 12 (Cn) Y n (Vn) 
2m^\Y 22 (( n ) Y 21 (r] n ) 



1 



det 



l e -f(V(Cn)-V(7?„)) 



x det 



/c a c« /2 ff2i(Cn) o N 



M(Cn) 



+ M( Vr , 



'0 r] a fjn 2 J a+ i(rjr< 
f'^J^ifc). 

"2 



We now rewrite the matrix in the determinant as was done in the proof of Theorem 11.11 Using 
also the fact that C a Cn H+i(Cn) is uniformly bounded for £ in compact subsets of C + as 

Q!±2 

n — > 00, see Lemma 14.31 we obtain instead of equation (|4.6I) . in a similar fashion, the following 

/rcn /2 ^ (1 ?i(c^ -^ 1/2 



1 



det 



V a T)n J a+ l{Vn) 



+ 0(l/n) 



(4.10) 



asn-> 00, uniformly for £ and r\ in compact subsets of C+ and C, respectively. We now insert 
the fact that, see Lemma 14.31 

CC l n /2 H^ ±1 _{L) = C a (fC) 1/2 ^iM + o(i/n), 



1 /2 

into (|4.1Uj) . and use the fact that rj~ a fj n J a +l{fj n ) is uniformly bounded for rj in compact subsets 
of C as n — > 00, see Lemma l4.2| to obtain 

^/,n +m (C^) = e-t^^)-^)) 

/C Q K) 1/2 ^iK) ^ Q ^ /2 J a+ i(r?„)^ 

r (l) , A s „-a~l/2 



idet 
4 



+ 0(l/n) 



(4.11) 



Inserting, see Lemma 14.21 

r ? -^y2j Q± i(r? n )=r / - Q (^) 1 /2j a± .(^)+0(l/n), 
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into P~TT|) . and using the fact that ^(ttQY^H^ 1 (jrC) ls uniformly bounded for £ in compact 

a± 2 

subsets of C+ as n — > oo, we then obtain 

Wn,n+m{(,V) 

U a H) l/2 H { ^Q r,- a (*v) 1/2 J a+ L(*r,? 



e 2 



1 , 

-det 
4 



+ 0(l/n) 



= e-S^)- v ^M ((C - vKjAC, V) + 0(1/")) , (4.12) 

as n -> oo, uniformly for £ and 77 in compact subsets of C+ and C, respectively. In (|4.12|) . 
J„ //(C) 7 ?) 1S l ne Bessel kernel given in Table Note that 

v'(o) ,f \ 

as n — > 00, 

uniformly for £ and 7/ in compact subsets of C+ and C, respectively. Inserting this into (|4.12|) 
and using the fact that (( — r])!^ //(C> v) ^ s bounded for ( and r\ in compact subsets of C + and 
C, respectively, the first part of the theorem is then proven. 

The second part of the theorem can be treated in the same way, using the behavior (|3.13|) 
of the second column of Y in the lower part of the disk Us, instead of in the upper part of the 
disk. □ 

We leave it as an exercise for the careful reader to prove the universal behavior at the origin 
of the spectrum for the kernel Wjjj, n +m- 

Proof of Theorem 11.31 The proof is similar to the proofs of Theorem 11.11 and Theorem 11.21 

□ 
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